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We investigate the orbital effects of strong external magnetic field on the ground state properties of 
a two dimensional holstein polaron, employing variational approaches based on the exact diagonal- 
ization (VAED). From the ground state energy and the wave function we calculate electron-phonon 
correlation function, the average phonon number and the Drude weight and investigate the evolution 
of a 2D holstein polaron as a function of the magnetic flux. Although the external magnetic field 
affects the polaron throughout the parameter regime, we show that the magnetic field has a stronger 
effect on a loosely bound (spatially extended) polaron. We also find that the magnetic field can be 
used as a tuning parameter, particularly for weakly coupled polaron, to reduce the spatial extent of 
a large polaron. 

PACS numbers: 74.50.-|-r, 74.20.Rp, 72.25.-b, 74.70.': 


I. INTRODUCTION 

The interplay between electronic and lattice degrees 
of freedom is central to many areas of condensed mat¬ 
ter physics. The holstein modeli, which is almost seven 
decades old, still holds a place of eminence, when it 
comes to electron-phonon (el-ph) interaction, because of 
the intricate many-body physics it encompasses within 
a simple paradigm. On the other hand, modification of 
electronic band structure on a lattice in the presence of 
strong magnetic field gives rise to the well-known Hof- 
stadter’s butterfly^, a quintessential example of fractals 
in quantum mechanical systems. Although the physics 
of solids is replete with electron-phonon coupled systems, 
precious little is known about the effects of strong mag¬ 
netic field on such a coupled system, where indeed only 
one of the partners, the electron, feels the effect of the 
field directly. 

In this paper we study the holstein polaron in a strong 
magnetic field, which marries the two issues above. In 
the presence of strong el-ph coupling, the electrons are 
expected to form polaronic bound states of varying size, 
depending on the strength of el-ph interaction. On the 
other hand, an orbital magnetic field induces a preces¬ 
sion of an electron away in an orbit reducing the net 
mobility, in effect subjecting it further to localization by 
phonons. It is in this context, the combined effect is an 
interesting puzzle calling for a resolution. However, the 
subject of lattice polarons in a strong magnetic field is 
largely an unexplored area, barring very few work, no¬ 
tably by Mona Berciii^, using momentum-average (MA) 
approximation. The variational approaches based on the 
exact diagonalization (VAED) used by Bonca, et &\^— 
and Chakraborty, et al4“— is one of the most successful 
numerical methods to study the holstein and extended- 
holstein type of el-ph systems in dilute regime in all di¬ 
mensions. 

We have generalized this method to deal with cases 
where the lattice unit cell has more than one equivalent 
site, i.e., a supercell-VAED. This scheme has been very 



FIG. 1: (Color online) A 2-dimensional square lattice in a 
magnetic field of The figure shows the variation of mag¬ 
netic flux along the y-bonds. 


successfully implemented in the present case and can be 
used to many important situations where a many-atom 
unit cell is coupled to another degree of freedom (like 
polarons in graphene or a case of partial disorder). We 
first compare our supercell-VAED zero field results to the 
benchmark results available in the literature, where we 
find at least a 8-digit match. Then we proceed to study 
the ground state Hofstadter band at different parameters 
of the holstein model. We first create a variational space 
by repeated action of the hamiltonian on the initial state 
and then we adopt a two fold approach: (i) integrate the 
spectral function obtained from k-space Green’s function 
over the brillouin-zone to get the density of states and 
(ii) find the ground state energy and wave-function by 
employing conjugate gradient technique^. 

This paper is organized as follows: in section II, we dis¬ 
cuss the hamiltonian and delineate the basis generation 
procedures. We then proceed to show our results and 
compare some of those with existing results in section 
III. In this section, we also study the evolution of the po¬ 
laron with magnetic flux in different parameter regimes 
and try to analyze the interplay between the el-ph in¬ 
teraction and the magnetic field. Conclusion follows in 
Section IV. 















2 


II. THE MODEL 

In order to study the effect of magnetic field, we employ 
the simplest and well-studied model for electron-phonon 
interaction, the holstein model, where a spinless electron 
is coupled to a dispersionless optical phonon, represented 
by a local Einstein oscillator. For a spinful model, the 
magnetic field leads to the Zeeman splitting of up and 
down-spin polaron bands. In a very high magnetic field 
and low electron density, one can consider only orbital 
effects and ignore the Zeeman splitting as the Hilbert 
space is spin-split into two sectors well separated in en¬ 
ergy. One can, therefore, concentrate on the lower of the 
two for low filling, effectively reducing the model to a 
spinless one. The hamiltonian is given by, 

H= -Y^tij{clcj + h.c)+ujY^blbi 

+ ( 1 ) 

i 

where < i,j > are near-neighbour site indices on a square 
lattice, Ci{bi) are electron (phonon) annhilation opera¬ 
tors respectively. The nearest neighbour hopping in¬ 
tegral in presence of magnetic field is now associated 
with a Peierls phase factor. The choice of a lan¬ 
dau gauge A{r) = B{0,ma,0) for a uniform mag¬ 
netic field B perpendicular to the plane of the lat¬ 
tice leads to the hopping integral tij = —t along x- 

direction and =—texp{±.2'Kim-^)= 

—t exp(±27rzTO|) along y-direction. Here (j)=Ba? is the 
number of flux quanta per plaquette, which is the gain of 
phase by an electron hopping round a closed path along 
the plaquette. ^ | with p, q co-prime integers and (po 

is the Dirac flux quantum. We set the hopping integral 
t to be 1 throughout the numerical calculation and all 
other parameters are defined in units of t. Here lo is the 
oscillator frequency and g is the dimensionless electron- 
phonon coupling strength. The effect of electron-phonon 
coupling is expressed in terms of two dimensionless pa¬ 
rameters a = LO It and g. 

Once the magnetic field is switched on, we lose lat¬ 
tice periodicity in x-direction. However, since it is only a 
phase, repeating at every 27r, periodicity is still retained, 
albeit, with a changed value. Fig.l, describes the situa¬ 
tion for a magnetic field B=^. In this case the changed 
periodicity of the lattice is 3, so instead of the one site 
unit cell, our magnetic unit cell is basically a 3-site strip. 
Therefore to account for a magnetic field B=^, the mag¬ 
netic supercell will be a strip of length N. The variational 
basis is similar to that of Bonca et. ah— , the difference 
being, now instead of one initial zero phonon state, we 
start with “N” zero-phonon states (for N different po¬ 
sitions of the magnetic strip) and while checking for the 
translational symmetry, we shift the supercell. Mathe¬ 
matically, this amounts to putting another index “m”, 




FIG. 2: (Color online) The upper panel (a) shows the lower 
lying landau levels for an electron in a lattice. The lower panel 
shows the same for a holstein polaron calculated at u;=0.5 and 
g=\/0.4:. The dashed vertical lines mark the ground state 
energy at Eq and Eq + u. 


to account for the position-dependent parameter. In or¬ 
der to get a convergent basis for small polaron (a;=5.0 
and g=2.0) we have appropriated the idea of Lang-firsov 
transformation while constructing the basis j^°’^^ff^ This 
method works very well for the case of small polarons as 
has been established by earlier works as it incorporates 
into the basis the important states which a small polaron 
requires 

We calculate the following quantities of interest: To 
calculate the correlation between the electron position 
and phonon distribution (lattice deformation) in the 
ground state, we define a correlation function as follows: 


x(*-j) = (V’gIcIg (fe](2) 

where ipG is the ground state wave-function. The total 
lattice deformation is conserved to 2g in the hamiltonian, 
from a straightforward sum-rule. The average phonon 
number is calculated as 
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(3) 


The distribution of number of phonons in the vicinity of 
the electron is given by 


7 (* - j) = ( V'gI c\ci IV'G ) (4) 

The Drude weight (-Do)of the ground state of the polaron 
is obtained by introducing a phase factor to the hopping 
matrix elements {t —> and then finding out the re¬ 

sponse to the electric current as 


Do = 


d^Eo{r])^ 
dr^2 l'^=o 


(5) 


where Eoif}) is the eigen energy of the ground state in 
presence of non zero rj^ The calculated Dq has been nor¬ 
malized with respect to a free electron on a square lattice 
for all the cases. 


III. RESULTS 

We first test the numerical accuracy of our result, by 
comparing the ground state energy of a holstein polaron 
at w = 2.0 and g = 1.0 with the best available results 
in the literatur o^d^ . The ground state energy for this 
parameter is —4.81473577, obtained from a vartiational 
basis constructed by operating the hamiltonian 11 times 
on the initial states (lV/i=ll— ), which matches up to 8 
decimal places with the result of Bonca et aX^. Then 
we compare the density of states (DOS) with the inclu¬ 
sion of el-ph coupling with the results of Berciu, et al4. 
Fig.2 displays the landau levels for a magnetic flux (^) 
of strength 0.005, for an electron (fig2.(a)) and that of a 
polaron (fig2.(b)) at w = 0.5 and g = \/0A (correspond¬ 
ing to A = 0.2 of Berciii^). In order to accommodate a 
flux of value 0.005, the size of our magnetic supercell (or 
magnetic strip)had to be 200. Our results are in excellent 
agreement with earlier results^ wherever a comparison is 
possible. The vertical dotted lines of fig 2(b) indicate the 
ground state energy {Eq) and Eq+oj respectively. The 
distinct polaron landau levels lose their identity beyond 
Eq+uj. We digress a bit and compare the scenario of 
an electron in a lattice with that of a holstein polaron 
without bringing the magnetic field into consideration. 
The electron has only a single band, whereas a holstein 
polaron has infinite number of bands as the electron is 
coupled to all possible phononic excitations. However, for 
holstein polaron we have a gap of w at k=0 in between 
the ground state and the first excited state, beyond that 
we have huge number of closely spaced states (depending 
on the parameter regime) and as we go up in energy they 
acquire a quasi-continuum nature. These are the states 


that disturb the sharpness of the landau levels beyond 
Eq+uj in fig.2(b). 

We now investigate the response of the polaron to mag¬ 
netic fields. The average phonon number (Nph) gives an 
idea about the phononic activity of the polaron. Increase 
in the value of Nph suggests an increase in phononic at¬ 
tributes of the polaron. Fig.3(a) shows the variation (rel¬ 
ative to zero-field value) of Nph at three different regimes. 
At w = 1.0 and g = 0.05, we have a quasi-free electron 
at i? = 0, and the phonon activity clearly increases with 
the increase in field values, achieving a maximum in the 
vicinity of </> = 0.375. Though the trend is similar at 
a;=2.0 and (/=1.00, the relative variation is much less 
pronounced. At w = 5.0 and g = 2.0, we are deep in the 
anti-adiabatic limit (w ^ t), polarons of small spatial 
extent hardly respond to the variation in magnetic field. 
Fig.3(b) shows the relative variation (w.r.t. the zero field 
value) of ground state energy with field for the same set 
of el — ph parameters. The pattern in ground state en¬ 
ergy variation is same for all the three parameter regimes 
and one can recognize them as ground state Hofstadter 
band. The relative change is again maximum for a weakly 
bound polaron and minimum for the polaron in the anti- 
adiabatic limit; the signature of magnetic field shows up 
in ground state energy pattern for all regimes. The varia¬ 
tion in average phonon number as well as ground state en¬ 
ergy with applied magnetic flux follow the lowest branch 
of the Hofstadter butterfly. Both fig.(3(a) and 3(b)) are 
symmetric about 0=0.5, a typical Hofstadter character¬ 
istics. The maximum size of the magnetic strip used for 
this calculation was 64. 



FIG. 3: (Color online) The upper panel shows the average 
phonon number Nph as a function of magnetic field for three 
different sets of el — ph coupling; Nph are normalized by their 
respective zero field values. The lower panel shows the ground 
state energy Eo (normalized by their respective zero field val¬ 
ues) as a function of magnetic field for three sets of el — ph 
coupling. 


A thorough study of the electron-lattice correlation 
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FIG. 4: (Color online) The electron-lattice correlation func¬ 
tion x{^ ~ y) rneasures the lattice distortion, shown at four 
different magnetic flux values. The el — ph parameter for the 
polaron is cu=1.0 and g=0.05. 


functions throws some light on the mechanism of in¬ 
creased phononic activity for a weakly bound polaron 
and a much smaller effect in the anti-adiabatic regime. 
A strip size of 16 was used to calculate this correlation 
function. Fig.4 shows the x(a;, y) for four different values 
of (j) {(f)=0.0, 0.1875,0.375 and 0.5) for a polaron at a;=1.0 
and 3=0.05. The displayed xi^tV) has been divided by 
the total deformation 2g for convenience. Though the 
total lattice distortion of the polaron always adds up to 
2g, we can see that with increase in y locally, the local 
distortion assumes a higher value for x(x, y) = x(0, 0) 
and its vicinity than farther outside, i.e., we a have a rel¬ 
atively tightly bound polaron (as compared to the 4>=^ 
polaron). The x(OiO) for (^ = 0 is 0.025 and reaches 
a maximum at </> = 0.375 to 0.043. A similar trend is 
observed for xix, y) as well. 

Fig.5 shows the electron-lattice correlation at same 
four values of (j) ((()=0.0, 0.1875, 0.375 and 0.5) for a po¬ 
laron in the anti-adiabatic regime (a;=5.0 and 3=2.0). An 
extremely small polaron results in this regime and as it 
is quite evident, almost the entire distortion is rooted at 
{x, y)={0, 0); there is no noticeable change with magnetic 
field. The x(O)O) for ()) = 0 is 3.90 and reaches a maxi¬ 
mum a.t (j) = 0.375 to 3.943, a very small change. A study 
of hg.4 and fig.5 clearly shows that the magnetic field 
brings about a prominent change in a weakly bound po¬ 
laron compared to a polaron that is tightly bound to the 
lattice. The magnetic field tends to shrink the distortion 
towards its centre for a weakly bound polaron, whereas it 
can hardly affect a strongly bound one. Clearly a loosely 
bound polaron has a much larger orbit and therefore has 
a stronger effect due to the orbital magnetic field. Hence, 
we see in fig.4, the relative increase in x(0,0), though the 
sum total of distortion was limited to 2g as expected. 

The Drude weight (do) gives the measure of coherence 
and is an important correlation function to study the na- 


FIG. 5: (Color online) The electron-lattice correlation func¬ 
tion x{^~y)y 8.t four different magnetic fields. The respective 
parameters are a;=5.0 and 3=2.0. 


ture of electronic conductivity of a system. There have 
been a number of earlier work suggesting magnetic field 
induced metal-insulator transition.— Fig.6 shows the cal¬ 
culated do as a function of field for the three el — ph 
regimes. The size of the magnetic strip used for this 
calculation was 16. At w=1.0 and 3=0.05 we have a 
quasi-free electron very weakly tied to the lattice, with 
a large spatial extent. Consequently, at <^=0, do ~ 1.0. 
However for very small magnetic flux (())=j^) it almost 
drops to 0, suggesting a complete loss of coherent hop¬ 
ping. The inset (a), (b) and (c) show the ground state 
band for five different flux values ((/)=0.0, ;|, | and i) 

for three different el-ph regimes. The inset (a), (b) and 
(c) at show a completely flat band (within the nu¬ 

merical accuracy of our calculation) and ^^=0 through¬ 
out the brillouin zone which arises from local quantum 
interferences— The electron, while moving through the 
lattice, will see a maximum variation in flux for small val¬ 
ues of field and loses its coherence, however the ground 
state band is not completely flat, for </> = ^ etc., 

we can find peaks in Dq. The Dq peaks for a;=2.0 and 
3=1.0, have slightly reduced values as we are in the inter¬ 
mediate el-ph coupling regime and we no longer have the 
quasi-free electron. For w=5.0 and 3=2.0, anti-adiabatic 
strong coupling regime, the exponential suppression of 
the do is seen, however the change in magnetic flux still 
has an effect, which can be seen in the inset (d) of fig. 
6 and local quantum interferences still has its effect at 
small magnetic flux {(j)=^), which leads to localization 
and consequent vanishing of Dq . We see a magnetic field 
induced coherent to incoherent transition at (j)= across 
the el-ph coupling regime. 































Drude Weight 
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FIG. 6: (Color online) Drude weight Do (normalized by the 
free electron Do) versus magnetic flux. The inset (a), (b) and 
(c) shows the ground state band at ())=0.0 (solid line), 0=]^ 
(dotted line), <i)=\ (open square), <)>=| (open diamond) and 
<()=i (dashed line) for (a;=1.0, <7=0.05), (a;=2.0, <7=1.0) and 
(a;=5.0, g=2.0) respectively. Inset (d) shows Do for ui = 5.0 
and <7 = 2.0 separately. 


IV. SUMMARY AND CONCLUSION 


We have developed a numerical scheme based on 
VAED to study el — ph interaction in presence of strong 
magnetic field. Our results are in excellent agreement 
with earlier results wherever a comparison is possible. 
The external magnetic field changes the orbit of motion 
of the electron and the effect is more in case of loosely 
bound polaron already spread out in real space. With 
the inclusion of magnetic field, we see that the spatial ex¬ 
tent of the polaron decreases resulting in a polaron with 
a tighter binding. Our method can be useful to many 
important cases with a more complex unit cell. 
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